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We prove on a prototype Hamiltonian that hidden long range order is always present in the gapped 
phases of interacting fermionic systems on one dimensional lattices. It is captured by correlation 
functions of appropriate nonlocal charge and/or spin operators, which remain asymptotically finite. 
The corresponding microscopic orders are classified. The results are confirmed by DMRG numerical 
simulation of the phase diagram of the extended Hubbard model, and of a Haldane insulator phase. 
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The behavior of strongly correlated electron systems 
has been widely investigated to understand the physics 
of several phenomena in condensed matter, ranging from 
the insulating regime to high-T c superconductivity. Due 
to the many degrees of freedom involved, many aspects 
of the micro- and macroscopic behavior of these systems 
remain unclear. Recently their simulation by means of ul- 
tracold gases of two-component fermionic atoms trapped 
onto optical lattices has opened new possibilities, lead- 
ing for instance to the direct observation of the predicted 
magnetic [T] and Mott insulating (MI) phases [2 . The 
latter is efficiently modeled by the Hubbard Hamiltonian. 
In this case, it has been noticed quite recently [3] that in 
one dimension (ID) it is possible to identify a nonlocal or- 
der parameter in the MI phase, which displays long-range 
order (LRO); a result that is in agreement with Coleman- 
Hohenberg-Mermin- Wagner theorem f¥| since no contin- 
uous symmetry of the system has been broken. The dis- 
covery envisaged a description of the underlying parity 
charge order, whose microscopic configurations are de- 
picted below in the second cartoon of Figjl] the Mott 
phase consists of a chain of single fermions with up and 
down spin, where fluctuations of pairs of empty and dou- 
bly occupied sites (holons and doublons) are bounded. 
The behavior is reminiscent of that observed in the insu- 
lating regime of a degenerate gas of bosonic atoms [5] . 

In general, the observation of gapped phases in ID 
systems is not believed to be necessarily related to the 
presence of some type of LRO, since the strong quan- 
tum fluctuation are expected to destroy any such order. 
In this Letter we show that LRO is instead hidden in ev- 
ery gapped phase of one dimensional correlated fermionic 
systems. The result is achieved by means of a general 
analysis of the bosonization treatment applied on a pro- 
totype lattice model Hamiltonian for these systems. We 
identify in the lattice the nonlocal parity and string op- 
erators responsible for the different types of LRO. As a 
byproduct, both charge and spin excitations turn out to 
be independently ordered, while local operators intrin- 
sically generate both. It is tempting to conclude that 
nonlocal operators are "more fundamental" with respect 
to the usual local ones, at least for the description of 



the possible orders in the ground state phase diagram of 
these systems. To test our results we perform a density 
matrix renormalization group (DMRG) analysis at half- 
filling and zero temperature of the standard extended 
Hubbard case, focusing on the insulating phases. 

We start from the general class of lattice model Hamil- 
tonians historically introduced in Ref.[fi] to describe the 
effects of Coulomb repulsion among electrons on their be- 
havior, the standard Hubbard model being the most fa- 
miliar example. The low energy behavior of these models 
is described by an effective Hamiltonian H obtained by 
bosonization treatment (see [7] and references therein). 
Upon neglecting terms of higher scaling dimension (see 
also [8]), T-L turns out to be the sum of two decoupled 
sine-Gordon models. Explicitly, we have 

n = E ( H o U) + 7^)2 [ dxco S [q v V8$ v (x)}\ , (1) 
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with ffW = %f dx[K„(jrTl u ) 2 + K- 1 ^,) 2 ]. Here 
is the compactified boson describing the charge {y = c) 
and spin (y = s) excitations, with velocity v v , Gaussian 
coupling K v and conjugate momentum n, y = OxQ^/tt: a 
is a cutoff. Moreover, in terms of the standard notation 
g c = g3±5 nrj -i 1 the corresponding term generated from 
Umklapp processes being non-vanishing only at commen- 
surate fillings n — p/q {p,q integer; we assume p = 1); 
g s = gix, and q c = q, q s =1. 

The cosine terms in Q become irrelevant in the renor- 
malization group (RG) flow equations unless the fields 
$„ are pinned to fixed values; in this case, the energy is 
minimized by the choices 

V2$ v = ^-(21 + 1) g„>0 (2) 

2q v 

V2®» = ~2Z 9, < (3) 

with I G NU{0}. A thorough inspection of the RG equa- 
tions shows that while both choices of locked values for 
$ c amount to the opening of a charge gap A c , a spin gap 
A s can open only for g s < 0, due to the SU(2) spin sym- 
metry of the Hubbard class of Hamiltonians. To resume, 
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Table I: Correspondence between ground state quantum 
phases and nonlocal operators that manifest LRO. 



in all systems described by T~L it is possible, depending on 
the filling factor and on the relevant interaction terms, 
to observe up to 6 phases (shown in Table |l| each one 
with different dominant correlations. In most phases the 
known dominant correlations of two-point local opera- 
tors decay to zero with distance following a power law, 
in agreement with bosonization predictions. On the one 
hand, in charge-density and bond-ordered wave (CDW 
and BOW respectively) phases - appearing when just on- 
site and nearest neighbors diagonal Coulomb interactions 
are present - LRO was identified with the non-vanishing 
in the asymptotic limit of appropriate two-point corre- 
lators of local operators [7 . On the other hand, quite 
recently it was noticed that for the standard Hubbard 
model LRO in MI and Luther Emery (LE) liquid phases 
is described instead by two-points correlators of suitable 
nonlocal operators [3]. In the present work we extend 
the idea of nonlocal order to all possible gapped phases 
listed in Table [T] for the general Hamiltonian H. 

First of all, we define for the lattice model the parity 
and string operators at a given site j as 
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vj„ is the number operator count- 



respectively, with v 
(rij^ — nj^). Here rij 
ing the electrons with spin a (a =t,i) at site j, namely 
rij a = Cj^Cje, Cj a being the operator which annihilates 

one electron of this type and cJ- CT its Hermitian conjugate; 
moreover rij = rij-f + rij^. The related two-point correla- 
tors Cp\r) = (Op\j)Op\j + r)) (parity correlator), 



and cf\r) = (0 [ g'(j)0 { g n (j + r)) (string correlator) 
can be approximated in the continuum limit according 
to the analysis outlined in Ref. OH], giving 

Cp (x) ps (cos V2$ u (0) cosV2$ 1/ (a;)) (5) 

iy/2* v (x)), (6) 



C { g\x) « (sm\/2$„(0) sin< 



where ( ) stands for the average evaluated in the ground 
state. From the above result one can realize that at least 
one of the parity or string correlators is non-vanishing for 
x — > oo in every gapped phase. Indeed these take place 
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Figure 1: Cartoon picture illustrating how order is maintained 
in presence of fluctuations. The blue continuos (dashed) lines 
show the correlated pairs of up-down spin (holon-doublon) 
allowing {0 ( p } ) ({0 ( p } )) to remain non vanishing. The green 
and red circles show the alternation of sites occupied by dou- 
blons and holons in the chain of single fermions preserving 
(Of) ft 0. 



when some <!>„ is pinned to a fixed values, as shown in 
Table|l] In that case, an order parameter (Oa^) emerges, 
since we observe 

lim C<M (x) = (OM > 2 = CW , a = P, S 



In Table |TJ LL stands for the gapless Luttinger Liquid 
phase, which is the only case without any LRO, as both 
the bosonic fields are unlocked (we indicate this with 
u) . LE is the conducting phase with open spin gap which 
takes place for $ s = 0, and is characterized by a nonzero 
(Op^). Charge-gapped phase with A s = can open for 
a) $ c = (MI), in which case (o£°) ^ [3 ; b) for 
c&c = 7r/v8, that we indicate as Haldane insulator (HI) 
since it is distinguished by the non-vanishing value of 
the Haldane-like string order (O^). Finally, BOW and 
CDW phases are fully gapped phases with two non van- 
ishing Oi^'s. Notice that only in these latter cases, the 
two nonlocal order parameters combine to form a true 
local LRO, namely the BOW and CDW orders that give 
the name to the corresponding phases [71 [S]. 

The non- vanishing of the parity and/or string correla- 
tors gives further physical insight about the kind of mi- 
croscopic orders underlying the phases. These are illus- 
trated schematically in Fig[T] At half-filling a non-zero 
value of the charge (spin) parity correlator implies the 
formation of bound pairs of holons and doublons (up and 
down spins) in a background of single electrons (holons 
and doublons) as it occurs in the MI (LE) phase [3]. 
Whereas a finite value of the charge (spin) string correla- 
tor, no matter the distribution of up and down electrons 
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(holons and doublons) , amounts to a holon (spin up) al- 
ways followed by a doublon (spin down) site on the holon- 
doublon (single electrons) sublattice. Not surprisingly, in 
this latter case the microscopic structure of the phase is 
quite resemblant of that of the HI observed in the bosonic 
case. Understanding the microscopic configurations in 
the different phases also allows insights onto the mech- 
anisms at the basis of the formation of charge and spin 
gaps. With respect to the perfect MI of singly occupied 
sites, the Mott charge gap at half-filling is maintained by 
adding localized pairs formed by a doublon and a holon; 
whereas the transition to an HI charge gap takes place 
when the added doublons and holons do organize into a 
sublattice obtained by ignoring singly occupied sites. The 
LE case illustrates how an open spin gap, ideally amount- 
ing to a configuration with holons and doublons only, is 
preserved when single electrons are arranged in localized 
pairs with up and down spins; this observation explains 
the reason for the dominant character of superconduct- 
ing correlations in such phase. Finally, combinations of 
the above possibilities determines the structures of the 
two fully gapped phases (CDW, and BOW). 

In order to support our predictions, we present below 
a numerical analysis of LRO parameters given by ^ , ^ 
for the insulating phases of the extended Hubbard model 
at half-filling in case of repulsive interactions. In this 
case the lattice Hamiltonian reads 
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Figure 2: Lower panel: Nonlocal order parameters C« (L/2) 
for L — > co for the three insulating phases of the extended 
Hubbard model Q. The vertical dashed lines locate the crit- 
ical points with the uncertainty determined from the numer- 
ical analysis reported in the upper panels. Upper-left panel: 
K s as a function of 1/L for different values of V. The con- 
tinuos lines are second order polynomial functions used to 
extrapolate the TDL. Upper-right panel: K c as a function of 
V for different system sizes. 



H = -*I>]> c i+i.* + H - c 

+ U n ft n 3l +V Y1 n 3 n i + 1 ( 7 ) 
3 3 

where U and V represent the diagonal on-site and neigh- 
bor sites contribution of the interaction potential and we 
fix the energy scale with t = 1. Such model, which is 
of fundamental relevance in condensed matter and in the 
younger field of ultracold systems, has been thoroughly 
explored in the literature (see [7] [TUJ QT] and references 
therein). Indeed, recent experiments with Fermi gas of 
magnetic atoms (T2] or polar molecules [13J allow to quan- 
titatively simulate the Hamiltonian ^ where both the 
interactions parameters can be tuned by changing the 
direction of the dipoles with external fields or the trans- 
verse frequency of the laser used to create the lattice. In 
particular, we want to explore at half-filling the regime of 
positive values of U and V , for which the phase diagram 
amounts to three different insulating phases. 
The analysis is performed using a DMRG algorithm on 
finite size chains with periodic boundary conditions. We 
have chosen to consider small system sizes, from L = 12 
to 48, with up to 1600 DMRG states and six sweeps in 
order to have a good precision on our quantities. 
We expect that the parity and string operators intro- 
duced above should work as order parameters for the 
three insulating phases. More precisely (see also Table IT|, 



the asymptotic value of (Op ) should be the only non- 
vanishing parameter for the MI phase; whereas in the 
BOW phase also (Op^) should become different from zero 
at the MI-BOW transition. Finally, at the BOW-CDW 
transition (O^) should become different from zero, while 
(Op) becomes vanishing. 

We have calculated C ( p ] (r) = (exp(£?=J mS ( " ] )) and 

C { s\r) = (SyeMYtt^l^S^sfl), and their 
asymptotic values have been evaluated at the mid point 
r = L/2, upon an extrapolation in the thermodynamic 
limit (TDL) L — > oo. Special care must be payed in 
separating the uniform and staggered parts of the parity 
operator, since the relation Cp\r) = (—l) r Cp\r) holds. 
Fig |2] collects our numerical results, showing a clear evi- 
dence of the expected behavior. 

Our findings can be compared with those obtained in fTl"] 
by considering the expectation value of a different non- 
local operator, namely the exponential position operator 
Zl- Since in bosonization analysis such value takes the 
form (cos V8$ c ), it is different from zero for both pinned 
values of <f> c allowed in an insulating phase, hence vanish- 
ing only at the conducting point where the BOW-CDW 
transition takes place |15j . 

To enforce our analysis we also computed the Luttinger 
constants K v defined as K v ~ lim g _>.o ttS u (q)/q, with 

S,(q) = iEu^ qik - l \(Sl z SL) - (SUWJ) in the 
TDL. These give precise information regarding the pres- 
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ence of gaps In particular the SDW-BOW belongs 
to the Berezinskii-Kosterlitz-Thouless (BTK) universal- 
ity class since a spin gap takes place entering in the fully 
gapped BOW phase, while maintaining a full rotational 
spin symmetry. The Luttinger theory predicts K s = 1 
in the gapless and K s = in the gapped phase. Nu- 
merically it is hard task to get exactly these two values 
since in the gapless phase logarithmic corrections affect 
the results, while in the gapped region really large sys- 
tem sizes are necessary in order to get K s = 0. It is 
customary to locate the transition point where K s takes 
values smaller than 1 in the TDL. As shown in Figj2] the 
transition point obtained in this way is in good agree- 
ment with the one predicted by O^ . The BOW-CDW 
transition requires particular care since its nature can be 
either second or first order, depending on the value of U. 
Here we consider the region U < 4 where the transition is 
known to be second order. As shown in [IT] , while the two 
phases are fully gapped, due to the competition between 
the onsite and nearest-neighbor interactions the charge 
gap is minimal at the transition point, where it takes the 
value 0. Hence the theory predicts a Luttinger parameter 
K c 7^ only at the gapless point and K c = elsewhere. 
In Fig(2]we show that K c develops a peak slightly depen- 
dent on the system size and extrapolations in the TDL 
confirm the transition obtained by looking the vanishing 
value of (Op?) and the finite value of (Og ). 

The scenario of Table IT] is completed by the identifi- 
cation of the HI phase, where only O s is predicted to 
have finite LRO. The ground state phase diagram of the 
model (|7| does not show such a phase [16_. Nevertheless, 
in Refsjll [17] a charge gapped phase corresponding to 
the pinned value $ c = was identified by adding 

to the Hamiltonian Q further correlated hopping terms 
of the form ^S(yy (T (cj .c,- (T + H.c.)(ra iff - n js ) 2 for an 
appropriate range of values of U and V. Such phase was 
denoted as bond-spin-density wave (BSDW), albeit the 
spin order cannot show LRO due to the unbroken SU(2) 
symmetry. On the basis of our analysis, since $ s is un- 
pinned, we expect such a phase to exhibit the searched 
HI order. We have numerically estimated the nonlocal 
correlators Ca (L/2) at various L in a single point in- 
side the phase (X = 0.25, U = 1, V = 0.5). The results 
shown in Fig|3] demonstrate that, within the numerical 
errors, in the asymptotic limit (and in the TDL) the only 

(c) 

operator that supports LRO is O s , as expected. 

We notice that further nonlocal orders may appear in 
fermionic systems as a consequence of reduced symme- 
tries. For instance, relaxing the SU(2) spin symmetry 
to U(l)xZ2, may allow for the appearance of the value 
$ s = 7r/\/8 in Eq.|l]), giving rise to Haldane-like corre- 
lations in the z-component of the spin. Further break- 
ing of the two U(l) symmetries related to particle num- 
ber conservation and spin rotation in the xy plane, open 
the way to a pinning of the dual fields O c and O s , re- 
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Figure 3: Nonlocal LRO in the Haldane insulator phase, at 
U = 1, V = 0.5 for the model ([7| plus a correlated hopping 
term with X = 0.25 (see text). As predicted from Table|TJall 
correlation functions Ca (L) vanish asymptotically except for 
Cg . Continuous lines represent nonlinear fits for estimating 
the asymptotic limit. 

spectively. As a consequence the correlators related to 
the operators cos(\/2©i>) and sin(v / 20 ! y) are also finite, 
thus generating a transverse Haldane-type order, simi- 
larly to what happens in spin-1 chains [18 or in the 
bosonic case [T51 I2U] . This simple argument suggests 
that for observing a Haldane order in all directions in 
fermionic systems it is imperative to extend interacting 
models like Eq.Q by including pair creation terms of the 
kind ^2j a (Cj a Cj +1 _ cr + H.c). Such analysis represents an 
intriguing topic per se that goes beyond the scopes of the 
present work and will be addressed elsewhere. 

In this Letter, we have proved that LRO underlies all 
the gapped phases of a large class of lattice model Hamil- 
tonians, describing ID correlated fermionic systems. The 
order is captured by appropriate two-points nonlocal cor- 
relators, describing separately the charge and spin excita- 
tions, which remain asymptotically finite in such phases. 
The observation allows to reconstruct the types of micro- 
scopic orders underlying the phases. In particular, it can 
be applied to discuss the phase diagram of the extended 
Hubbard model at half-filling, as well as the presence of 
an Haldane insulating phase in fermionic systems. Our 
results give precise indications as for the nonlocal quan- 
tities to detect in experiments with fermionic trapped 
dipolar atoms. These are directly accessible to experi- 
mental detection in optical lattices via single site resolu- 
tion imaging with the available techniques [5]. 

The generality of the analysis here described suggests 
the presence of a universal mechanism extendable to any 
system in ID, stating the presence of appropriate nonlo- 
cal LRO in every phase that shows a gap in the excita- 
tion spectrum. A further interesting topic which is still 
under debate concerns the relationship between the role 
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of nonlocal order, topological phases, and long distance 
entanglement [5T] in ID. 

Regarding the possible presence of the discussed type 
of nonlocal orders in higher dimension, one should stress 
the different nature of the parity versus string order pa- 
rameters, as depicted also in the cartoons of Fig[T] In- 
deed, the parity LRO is preserved by adding pairs of par- 
ticles that locally break the order, and this can be done 
in arbitrary dimension. At variance, the string order re- 
quires that the added particles organise into a nonlocal 
ID structure, a feature that seems difficult to generalise 
to higher dimension. With this in mind, it is reasonable 
to expect that phases with just parity order parameters 
(MI, LE, and BOW) could describe the physics of two- 
dimensional systems as well. This is consistent with very 
recent results on the relevance of parity correlator in the 
MI phase of the 2D Bose-Hubbard model [22 , as well as 
with findings on backflow correlations in the 2D Hubbard 
model [23 , which emphasize the role of holon-doublon 
attraction in the MI phase. 
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